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Abstract 

We examine in non-Abelian gauge theory the heavy quark limit in the presence of 
the (anti-)self-dual homogeneous background field and see that a confining potential 
emerges, consistent with the Wilson criterion, although the potential is quadratic and 
not linear in the quark separation. This builds upon the well-known feature that propa- 
gators in such a background field are entire functions. The way in which deconfinement 
can occur at finite temperature is then studied in the static temporal gauge by calcula- 
tion of the effective potential at high temperature. Finally we discuss the problems to 
be surmounted in setting up the calculation of the effective potential nonperturbatively 
on the lattice. 
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1 Introduction 



Over the years, various characterisations have been proposed for 'confinement', the prop- 
erty that coloured degrees of freedom are undetectable at present-day collider energies. 
Certainly, the Wilson criterion that static colour sources cannot be separated arbitrarily 
far apart [||] has lead to many insights both in lattice simulations and in analytical calcu- 
lations. In particular, based on the Wilson criterion lattice simulations have established a 
confinement-deconfinement phase transition at finite temperature. There are however alter- 
nate characterisations for confinement which may be more directly relevant for dynamical 
quarks and gluons, and which are based on the analytic properties of the nonperturbative 
quark or gluon propagators. In this paper, we shall focus on the suggestion that the absence 
of poles in the complex energy plane of field propagators is consistent with confinement of 
quarks and gluons, in other words that propagators are entire functions. That this can be 
correctly described as 'confinement' is easy to see: absence of poles means that no coloured 
degrees of freedom can appear in physical asymptotic states. This characterisation of con- 
finement is not necessarily in conflict with the Wilson criterion. Indeed, one of our aims 
will be to show that, in the static quark limit, entire quark propagators lead to the Wilson 
criterion. 

A quite simple mechanism for rendering quark and gluon propagators entire in the 
complex energy plane is to apply a homogeneous background gluon field which satisfies the 
key property that it be either self-dual or anti-self-dual. Such a background gauge field is 
characterised by 

5^(x)t'^ ^ ^nV^S^^x, (1) 

Bfj.u = -^e^upxBpx = ±Bp^ (2) 
B^uBpp = B'^Sup, B = const, 

Bij = —CijkBk, Bj4 = ±Bj . (3) 

The positive and negative sign in Eqs.(0j3|) correspond, respectively, to the self-dual and 
anti-self-dual cases. The colour vector points in some fixed direction which can be 
chosen such that n^r"^ is diagonal; n"" picks out the Cartan subalgebra of the colour group. 
Various properties of this field in SU(2) gauge theory were investigated originally in [§, ^. 
For example, in contradistinction to the chromomagnetic background field 0| the self-dual 
background is stable. Moreover, it was observed that this field leads to entire functions for 
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the charged scalar field propagator. In the sense described above then, this field can provide 
for confinement of quarks and gluons. Diagonal components of the gluon field (as SU(N) 
algebra elements) are not confined at least at the level of the lowest order propagator in 
the background field. A self-dual homogeneous field is at least then a possible source for 
confinement in QCD if it can be shown that such a field is a dominant configuration in the 
QCD functional integral. 

This verification can come from a computation of the effective potential for the candidate 
background field and the demonstration that the potential has a minimum at a nonzero value 
for the background field. The effective potential was calculated to one-loop in |§, ^. These 
results however were inconclusive in the sense that the quantum corrections to the potential 
were as large as the zeroeth order classical term. To our knowledge, despite several attempts 
to study the effective potential for the Savvidy chromomagnetic background |Q on the lattice 
lli 01 li 0> IC], an analogous nonperturbative computation for the self-dual homogeneous 



background has not been attempted. Nonetheless, with the assumption that the effective 
potential for a self-dual homogeneous background field has a nontrivial minimum and using 
just those quark and gluon propagators which exhibit confinement in the sense of entire 
functions, some successful phenomenological investigations for SU(3) have been carried out 
|11, Quantitatively, experimental data for the spectrum of light, heavy-light and heavy 
quarkonium systems can be reproduced to within ten percent in this effective description. 

In this work we concentrate on the problem of confinement and the effective potential 
for the SU(2) gauge theory. Our goal is first to describe the confining properties of the 
self-dual background field in the more familiar terms of the Wilson picture Secondly, 
we seek to show that, even if we cannot prove the existence of a nontrivial minimum in 
the effective potential for this background field at zero temperature and strong coupling, 
nonetheless deconfinement at high temperature can occur. Namely, we will show that at high 
temperature the effective potential for a self-dual background field acquires a minimum at 
zero field value. 

In the first instance we illustrate the confining properties of the self-dual homogeneous 
background by studying the problem of heavy particle and anti-particle in this background 
field. We thus examine the nonrelativistic limit. We indeed find that a confining potential 
for static charges emerges: the stationary trajectories of particle and anti-particle in the 
background field (||) separated by distance \X\ and held apart for time T are suppressed by 
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a factor 

exp 

where is a reduced mass of two-particle system. This result differs from that seen in most 
lattice simulations because of the different long-range properties of the field considered here 
as compared to those normally implemented in lattice gauge theory. The oscillator binding 
potential arises here effectively due to an interaction of the charges with the background 
field, but not by virtue of quantum gluon exchange between these charges. The self-duality 
and homogeneity of the background field is of crucial importance. The oscillator nonrela- 
tivistic potential is not inconsistent with the phenomenology of Regge trajectories in the 
hadronic spectrum since the latter is a feature of light quark systems. In the approach 



to the relativistic bound state problem of |12] based on the bosonisation of the one-gluon 
exchange interaction between quark currents in the presence of the vacuum field @) it is 
seen that the property that quark and gluon propagators be entire precisely gives rise to 
Regge behaviour in light-quark systems. 

It is appropriate to mention here the evident fact that a vacuum field such as Eq.(|l|) 
would lead to breaking of a range of symmetries such as CP, colour and 0(3). A satisfactory 
restoration of these symmetries at the hadronic scale assumes the inclusion of domain struc- 
tures in the vacuum. In a given domain the vacuum field has a specific direction and is either 
self-dual or anti-self-dual, but this is uncorrelated with the specific realisation of Eq.(||) in 
another domain, configuration The idea of domains in the QCD vacuum was discussed in 



application to various homogeneous fields IQ, 13, In the effective meson Lagrangian 
of |jl2| this idea was realized as the prescription that different quark loops (namely, those 
separated by the meson lines) in a diagram must be averaged over different configurations 
of the vacuum field (|l]) independently of each other. In the present paper we do not con- 
sider this problem, and only wish to note that the above formula for the contribution of 
stationary trajectories does not depend on directions and is the same for both self-dual and 
anti-self-dual homogeneous fields. 

In the second instance, though we cannot compute the effective potential nonpertur- 
batively we nonetheless seek to show that at high temperature, where asymptotic freedom 
should set in, the effective potential does actually acquire a minimum at zero external field 
consistent with deconfinement. This is not just a trivial consequence of perturbation theory. 
Lattice simulations have confirmed the picture that high temperature Yang-Mills theory. 



though deconfined, shows significant signals of nonperturbative structure In order to 
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account for some of these properties we have used the recent developments in temporal/axial 
type gauges at finite extension or temperature by Lenz et al fl^, 17]. Here a complete gauge 



fixing of Yang-Mills theory was formulated, accompanied by an integration out of certain 
zero mode fields which themselves are related intimately to the Polyakov loop order param- 



eter 1 18] for the confinement-deconfinement phase transition in pure Yang-Mills theory. The 
integration out of these variables generates for off-diagonal gluon fields a temperature de- 
pendent mass M{T) which diverges with increasing temperature, T. In ]|l^] it was checked 
that, despite the gluon mass, renormalisation at the one-loop order was standard, leading 
to the correct one-loop beta function for SU(2) consistent with gauge invariance. Moreover, 
this mass was shown to be related to the string constant in a linearly confining potential. 
Though the actual mechanism for confinement in our study is quite independent of that in 



]16], this gluon mass generation is of crucial importance for us. It defines a scale M = M{T) 
in the running coupling constant g^{M) so that at high temperature the coupling is small. 
We are thus able to perform a controlled calculation and find that at high temperature the 
effective potential takes the form, 

U^B^) = ^7?^ + o f X + OiB^/M^iT)) + OigliM)) 
^ ' 5|(M) 5257r2 M4(T) ^ ' ^ " ^^^^ " 

which has a minimum at zero field 5 = 0. If non-zero B can generate confinement at zero 
and low temperatures, then our result shows that deconfinement at high temperature can 
occur. 

In the following section we demonstrate that the self-dual homogeneous field provides 
simultaneously for the Wilson confinement criterion and the property that propagators of 
off-diagonal (charged) fields in a self-dual homogeneous gauge field are entire functions. 
Following that we consider the high temperature limit in the effective potential. The paper 
concludes with a summary of results and a discussion of the problem of computing the 
effective potential on the lattice. Much of the detail of explicit calculations is relegated to 
four Appendices. 

2 Self-dual Homogeneous Field and the Wilson Criterion 

To illustrate the relationship between confinement and the property that Green's functions 
in an (anti-)self-dual background field are entire functions it suffices to consider a simple 
charged scalar field of mass m coupled to the background gauge field = B^^Xy defined 
by Eqs.(||-0). The relationship between this and the original Yang-Mills theory can be 
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understood as follows: by assumption, the effective potential for the configuration Eqs.(|^ 
|3|) exhibits a minimum at ^ which itself is proportional to the fundamental scale of 
the theory, Aym- By shifting the fields, we study the coupling of small fluctuations to this 
non- vanishing background. Thus the (/)-fields are those components of the gluon field which 
couple in the leading order to the background. We are thus lead to the effective Lagrangian 

and work, initially at least, in Euclidean space. Because we seek to approach the Wilson 
criterion, we consider the analogous Green's function describing a particle-antiparticle loop. 
Thus the object we are interested in is the four-point function 



G{x,y\B) = (: cp^x)(^ix) :: <t>HyMy) :)b = S{x,y\B)S{y,x\B) 



(4) 



The normal ordering is taken to exclude the disconnected diagram. The two-point function 
S{x, y\B) is itself a solution to the equation 

'-{d^ + iB^{x)f + m^] S{x,y\B)=5{x-y). 

The propagator in the external field transforms under translations {x ^ x + a, y ^ y -\- a) 



as 



S{x, y\B) = et^''-^'^''"''5(x + a,y + a|5)e~i^^^''-''-. 



(5) 



The Green's function Eq.(^ is gauge invariant and, hence, translation invariant. By means 
of transformation Eq.(|5|) with a = —{x + y) /2 we rewrite the function Eq.(^) in a manifestly 
translation invariant form: 



G(x, y\B) = G{x + a,y + a\B) = G {{x - y)/2, {y - x)/2\B) = Wix - y\B). 

Using the proper time method the propagator can be represented in the form of a path 
integral over a one-dimensional field ^ JH 



S{x,y\B) = e 



, j DeexpJ "I dr^ i\T)+ii^{T)B^,Wr) ) (6) 











with the boundary conditions ^(0) = — (x — y)/2, ^(a) = {x — y)/2, and the normalisation 
j Deexpj- Jdr^l =exp{-(x-y)V2«}. 
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Let us first review the confining properties of these fields in terms of analytical properties of 
the propagator. It is instructive to consider first the case of arbitrary constant B^,^. Since 
vectors H ± E (Hi = E-i = B^) are rotated independently of each other under 

Euclidean 0(4) transformations, the tensor B^u can be put into the configuration B34 = E, 
B12 = H, Bi3 = Bu = B23 = B24 = 0, and H > 0, -H < E < H ^]. The path integral in 
Eq.(^ can be easily performed with the result 

00 9 
H\E\ f dae-"" " 







sinh{Ha) sinh(|i?|a) 
X exp l^-^H[{xi - yif + (xa - y2f] coth{Ha) 

-^1^1 [(X3 - ysf + - y^f] coth(|E|a)| . 

This leads to a Fourier transform of the translation invariant part 

00 9 

"^^^'^^ " y cosh(i/a)cosh(|^|a) 

X exp [-j^{pI +P2)tanh(i/a) - ^(^i tanh(|^|a)| . (7) 

When E is nonzero this function is finite for any complex p\ + ^-nd p\ + p\ and thus is an 
entire analytical function. When E = Q this representation exhibits a pole in the physical 
region p\ = —{p\ + m? + H), which corresponds to a free propagation along the third axis 
with the energy equal to the lowest Landau level of spinless particle. In the (1 — 2) plane 
the particle is confined. 

Thus for 7^ 0, no physical particle corresponding to the field (t){x) can appear in 
the spectrum. The charged particles are, in other words, confined. However, as has been 
shown in |^], such an abelian constant field is unstable against small quantum fluctuations 
except in the case that it is self-dual or antis elf- dual: H = B, E = zizB. In the following, 
we concentrate precisely on this configuration. In this case Eq.(0) takes the simple form 
(i = tanh(5a)) 

S{P'\B) = ^ dtr—1] exp -^t (8) 



BJ \l + t, 

^ ' 

which represents an entire function in the complex p'^ plane. A special case is that of m = 0: 
the Fourier transform of the massless propagator turns out to be 

S{p'\B)\m=o = (l-e-P'/'')/p' . (9) 
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This is manifestly an entire analytical function in the complex p^-plane: the apparent mass- 
less pole at = simply cancels out illustrating most cleanly the confinement property. 
As a matter of fact, entire propagators mean that the quantum field theory is nonlocal. It 
should be noted here that at the axiomatic level nonlocal quantum field theory was suc- 



cessfully constructed some time ago [EQ, BT^, 22|. In particular, causality and unitarity of 



the S"— matrix was proved, a procedure for canonical quantization of nonlocal field theories 
was constructed and, recently, Froissart type bounds on cross-sections at high energy were 



obtained |23|. But to summarise this brief review of known results for constant fields, we 
can say that confinement in the sense of entire propagators is a property of any Euclidean 
abelian constant field configuration with non-zero magnetic and electric components, but 
the (anti-)self-dual case is distinguished by being stable against quantum fiuctuations. 

To see how this property can relate to the Wilson criterion, we now approach the problem 
of static charges. We consider heavy particles, with ^ B. In this limit Eq.(P) can be 
represented in the form of a quantum mechanical path integral (see Appendix A) 

S{x, y\B) oc e-'"^ J Dffexp | - ^ (i/3L(r/(/?)) | , (10) 

where 

T = X4-y4, f/(0) = -(x - y)/2, ff{T) = {x - y)/2. 

Here, Ej = B^j is the electric component of the tensor -B^i/, and Bi = —\e-ijkBjk is the 
magnetic component. We will implement the (anti-)self-duality condition Ej = iBj below. 
For the present, we insert the representation Eq.(^) into Eq.(^, introduce the center of 
mass coordinates R = {ffi + f]2)/2, f = fji — f]2, R{0) = R{T) = 0, r(0) = —riT) = y — x, 
and integrate out the center of mass coordinate R. The integral over R obviously does 
not depend on x and y, which is simply a consequence of the translation invariance of the 
function W. After continuation to physical time (T = iT, [3 = it) the result for W is 

W{x-y,T\B) oc e-^'""^ J Dfexpl^i J dtL{f{t)) 

■-^ -I -I 

^ = ^ + -/[r>,B\--{rB)\ 
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where = m/2 is the reduced mass of the two-particle system. One sees that the conjugate 
momentum and the Hamiltonian are 

and that the function W can be reexpressed as a phase-space functional integral, 

W{x-y,T\B) oc e"^'""^ J DrDp exp |-i J dt \H{r,p)-p- f] | . (12) 

Eqs. (|l^) and (|l^) show that the massive charged particle and anti-particle in the external 
self-dual field are bounded by an oscillator potential. Now, consistent with Wilson [Q], 
we extract from the path integral the contribution to the phase space of the stationary 
trajectory (p = 0, |f| = |x — ?/|). Equation ( [Tl| ) indicates that this trajectory corresponds to 
uniform circular movement of the particle-antiparticle pair on a circle with radius \x — y\ 
in the plane perpendicular to the direction of field B. We find that the contribution is 
exponentially suppressed, 

exp|^-ir|^(x-y)2^ . (13) 

The Wilson criterion for confinement is indeed satisfied. However here we have a 'volume 
law' rather than an area law. The relationship between this result and that in standard 
lattice gauge theory will be discussed in the final section. For now, we stress that the 
confining potential has appeared due to the background field, and not due to interaction 
between particles via gauge boson exchange. Such effects will generate additional potential 
terms to the Hamiltonian, and will thus affect the energy spectrum of the system. But gauge 
boson exchange will not change the basic confining properties of the background field. 
This picture of bound state formation seems strange at first sight. However, an analogy 



with the quantum dots (or artificial atoms) of solid state physics can be recognised |24]. 
Quantum dots are quasi-zero-dimensional electron systems in semiconductor nanostructures 
in which three-dimensional confinement of small numbers of electrons is achieved by a 
combination of band offsets and electrostatic means. The simplest model Hamiltonian 
for the few-electron quantum dot was obtained by solving the Schrodinger and Poisson 
equations self-consistently within the Hartree approximation [^]. It was found that the 
oscillator confining potential has nearly circular symmetry. The difference in our case is 
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the origin of the confining potential. The Hamiltonian Eq.(ll) has appeared due to the 
background gauge field which may arise in the vacuum as a result of gluon self-interactions. 
In QCD, this picture of confinement and bound state formation in the static quark limit 



will be basically the same. Thus Eqs.(ll,12) give illustrative insight into the basic nature of 
confinement provided for by the self-dual field. But, as Eq.(^ indicates, the significance of 
the property of entireness of Green's functions as a characterisation of confinement applies 
to dynamical fields and thus is relevant to the fully relativistic bound state spectrum of 
QCD, the physically relevant problem. Thus the qualitative basis for investigation into the 
impact of confinement on the relativistic bound state spectrum are equations like Eqs. (jsl) 



and (y), as has been carried out in |12]. Here an effective meson theory based on the 
bosonisation of nonlocal quark currents has been developed. The background field has been 
taken into account both in quark and gluon propagators. Within this effective theory the 
ground and excited state spectra of light, heavy-light mesons and heavy quarkonia have been 
calculated, with the only parameters being quark masses, the background field strength and 
the gauge coupling constant. Agreement with experimental data is obtained to within ten 
percent. Regge behaviour within this approach is recovered precisely by the fact that gluon 
and quark propagators are entire functions. The relationship between this mechanism of 



confinement and flavour chiral symmetry breaking is analysed in |11 



Having explored again the confining properties of the self-dual homogeneous background 
field in QCD, we now turn to the problem of the effective potential for this field at finite 
temperature and the question of deconfinement. 



3 Self-Dual Field and Finite Temperature 

In this section we compute the one-loop effective potential for the self-dual background field 
at finite temperature in SU(2) Yang-Mills theory. This enables us to study its presence or 
absence at high temperatures where perturbation theory should become reliable. 

Since we are already in Euclidean space in order to define the self-dual field, it is conve- 
nient to introduce finite temperature T by working in the imaginary time formalism. The 
X4 direction is now a finite interval of length (3 = 1/T and boundary conditions must be 
imposed on the gluon fields to which we shall come below. We work in a completely gauge- 
fixed formalism within which we will introduce the external field. At zero temperature, the 
background gauge is most convenient. However in the present case, the breaking of manifest 
Lorentz invariance (by the heat bath) suggests the temporal (axial) gauge is a natural gauge 
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choice. Specifically, we choose 

d^AKx) = (14) 

followed by a diagonalisation of the surviving zero mode a2{x)T"' = ^ Jq A2{x)T°'dx4. This 
gauge is a special case of the static temporal gauge. Here one encounters the problem 
of the nontrivial Haar measure in the functional integral quantisation of the theory p7[ |. 
Concomitantly, the diagonalised variable a^^'^^{x) is compact. The functional integral over 
this variable is thus non-Gaussian. Progress on the computation of this integral for SU(2) 
Yang-Mills theory was made recently in |jl^ wherein, using a lattice regularisation, it was 
shown that the integration out of a'l^^^{x) leads to an effective action for the remaining 
degrees of freedom. In the absence of external fields, the key features of this effective 
theory were that off-diagonal, namely charged, components of the gluon fields acquired a 
temperature dependent mass M(T). Secondly, the boundary conditions in X4 of these fields 
were changed from periodic to antiperiodic. 

We rederive this effective theory in Appendix B, and show that the presence of the self- 
dual background field does not force major modifications. In particular, the rigorous result 



for the mass, expected to be valid at low but non-zero temperatures [16, 17] is reproduced 
even in the presence of the homogeneous field, namely. 



M(T) = J (7r2/3- 2)r . (15) 



In [17| it was argued that stability with respect to chromomagnetic fluctuations mean that 



the mass term in the deconfined phase should take the form 

M{T) = ^Tg^{T), T^oc , (16) 

where g{T) is the perturbative running coupling constant. The important consequence 
of this result is that at high temperature the mass itself diverges but the ratio M{T)/T 
vanishes in this limit. This latter property is sufficient to guarantee the recovery of the 
Stefan-Boltzmann law in the high-temperature regime. 

Now we consider the self-dual external field and choose it to point in the same colour 
direction as air"". It is important to note that this corresponds to a distinct physical choice 
since gauge freedom does not allow both i?^r" and a^r'^ to be simultaneously diagonal. 

We come to the question of the gluonic boundary conditions. Here care is required 



as, unlike the chromomagnetic choice Q 26, 17], the self-dual field involves a component 



pointing in the, now compact, time direction. We are therefore no longer free to impose 
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the usual periodic boundary condition. Instead, the choice must be consistent now with 
parallel transport in the presence of an external field. Specifically, the appropriate boundary 
condition in the spatial directions x is the usual vanishing one. For the direction X4 which 
is finite, X4 S [0,(3], one usually chooses periodic boundary conditions in the absence of 
external fields. This can be represented in the form 

e^^^A;i(x4,x)=^;i(x4,f). (17) 

In the presence of an external field the natural generalisation of this for the fluctuating 
gauge fields is obtained via parallel transport, namely 

{e^'''y''Q''^ix,,x) = Ql{x^,x) 

Df = 6"'' 84 - e^^^B^ . (18) 

This boundary condition will, in the simplest way, preserve the periodicity of observable, 
gauge invariant quantities. We shall refer to this position-dependent twisted boundary 
condition as quasiperiodic. When the considerations of Appendix B are carried out and 
the zero mode of the fluctuating gauge field, q^^^^, is integrated out, Eq.([l^) becomes a 
quasi- aniiperiodic boundary condition: 

(e^^*)"''Q^(x4,f) = - Q;1(x4,x) . (19) 

To summarise what will be important then for the following calculation, there are two 
key features: firstly, that boundary conditions are modified to being quasi-antiperiodic, 
and secondly that the off-diagonal gluon components have a temperature dependent mass 
M{T) which diverges as T increases. It is precisely this which gives us a well-controlled 
high temperature regime specified by T S> Agu(2) and B < T^. 

To calculate the effective potential now, it is convenient to bring the field-strength tensor 
to the form (taking the field B to be directed along the third spatial axis) 



^l,ly=l,2,3,4 



/ -B \ 

BO 

±B 

V ^B J 



(20) 



where the upper (lower) sign corresponds to the self-dual (anti-self-dual) field. The effective 
potential is defined in the usual way using the functional integral. 



J DQexpl^J d^xCes[Qt,QlBf,]^ = exp{-(3VU,s{B,P,g)} (21) 
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CcE[Qt,QlBf,] = £ym[Qm>^m]Iq4=o - ^M\T)QtQ 



A 
i 1 



(22) 



where i,j,k,l = 1,2,3, A,B = 1,2 denote spatial and off-diagonal field components for 
gluons respectively, V is the three-dimensional spatial volume, and, as derived in Appendix 
B, the effective Lagrangian can be written as 

1 

2' 

with £ym the standard Yang-Mills action. The functional integral is defined on the space 
of quasi-antiperiodic fields satisfying Eq.(|l9|). The normalisation in Eq.(pTD is chosen so 
that Uef[{0, P, g) = 0. To the action, a gauge-fixing term involving the neutral zero mode 
gluons Qi{x) = ^ Jq dx4Qf{x) can be added, but which decouples at one-loop after the 
normalisation at zero field, B = 0. Dropping terms in the Lagrangian higher than quadratic 
in the fluctuating flelds Qf, we can extract from Cym the following piece relevant for the 
one-loop effective potential 



C 



1 



Qt{x) -{vY^'^ij + M^(r)5^^5,, + {D,Dj 



AB ! 



\AB 



+ 2Bije 



3AB 



Qf{x) (23) 



where = D^ + Df. The quadratic operator in Eq.(^) has zero modes for the case M = 
which are called chromons 0. A correct calculation of the chromon contribution to the 
effective potential at zero temperature requires an extension of the one-loop approxima- 
tion: an interaction (or mixing) between zero modes and normal modes has to be taken 
into account. However, we take the temperature to be sufficiently large so that M^(r) is 
correspondingly large compared to B. Large M means that the contribution of chromons is 
regular at one-loop order so the mixing between them and normal modes can be neglected. 
The one-loop effective potential is thus given by 



B^ 



1 



Ucs = ^ + TTT^Trln 



{-V^Sij + DiDj - 2iBij + M^{T)5ij) 



(24) 



gl ' V(3""' i-d^ki+dkdi + M^T)6kU 

Here B^/gQ comes from the classical action with go is the bare coupling constant. The 
effective potential can be rewritten in the form 



B 

Uea = — 2" 
9o 



(fx 



dX4 



dm 



D^Ax, x\B, m) — D'/Ax, x\0, m 



Pi 



(25) 



We have to calculate the trace of the propagator D^-{x, y\B,m) satisfying quasi-antiperiodic 



boundary conditions, Eq.(|l9D: 

D'^-{xi + (5,x]y\B,m) 
DiAx;y4 + P,y\B,m) 



-D^Axi, x;y\B,m) exp[-f/3i?4(x)]. 



D'-Ax; y4,y\B,m) exp [2/3^4 (y)]. 



(26) 
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This can be implemented by first solving for the Green's function Aij{x\B,m) relevant to 
the zero temperature or infinite volume and then building up the Green's function satisfying 
the finite temperature boundary condition via (see also [E^] and references therein) 



Dij{x;y\B,m) 



{-!)"■ Aij{xi-yi + nP;x-y\B,m) 

n=—OD 

X exp (^^Xf,Bf,iyy^ + ^nPBi{x + y)^ . 



(27) 



It should be stressed that Eq. ( ^7|) implies the existence of an orthogonal complete set of 
eigenfunctions of the operator V'^ satisfying the quasi-antiperiodic boundary conditions. 
The existence of such a set of functions is demonstrated in Appendix C. 

The infinite volume or zero temperature Green's function Aij is a solution to the equa- 
tion 



+ I 6ij + DiD 



2ifi 



Ajk{x\B,m) = 5ik5{x) 



(28) 



A complete solution of this system of equations is quite involved, but the trace of the 
propagator is tractable as is shown in Appendix D. One comment is in order though: the 
summation over n in the space-time trace of Eq.(p7|) is suppressed in the infinite volume 
limit V ^ CO due to the electric field component of the self-dual field. So in fact the only 



relevant contribution of Eq.(|27|) to the effective potential Eq.(25) is that from n = 0. Using 
this fact and results Eqs. (p.5|) and (p.2|) derived in Appendix D, we arrive at the relation 



An(0|B,m) 



E 

k 



Fk{x,x) 



(fzFk{z,x)bl{x)Ai{x,z) 



+2B^ J (fz J d^z'Fo{x,z)Ai{z,z')Fo{z',x), 



(29) 



where 

Fk(.x,y) 

A4{x,z) 



fi ^ \ B^ dr 



X exp 
1 



(Br) 



-m^r + 2B^kr - -{x - yfBcoth{Br 



(30) 



2^^ 



dt 

71 



exp 



2 {Z4 - X4,y 

-m t — 



At 



{Z4 — X4)BijX^ 



According to Eqs. (p5[), ( p9| ) and (|30|), the effective potential can be expressed as the 
combination 



f/cfr(i?' 



B' 



9l{M) 



+ Ui{B'') + U2{B^) + U^{B^), 



(31) 
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where 




^2 



^1 = / — exp( — — g ) I . [1 + 2 cosh(2s)] - 3s' 

R2 I- i-oo rlarli ( 

U2 = 




+ 



327r2 77o s2(s + t) 

f 2 smh(2s) - coth s [1 + 2 cosh(2s)] 

I sinh^ s \/l +t coth s 

3 ts2 ^ 



X 



I [sinh(s + r)] [sinh(s + r) + i cosh(s — r)] 



-1/2 



s + r 



-3/2 



(s + r + t)-i/2}. (32) 



The functions C/i, U2 and f/a correspond to ultraviolet finite contributions of the three terms 
in Eq. (]29|). The renormalised coupling constant g]^ runs with the scale defined by the mass 
M = M{T) and is 

4 (i_,„r (33) 

9r % \ Jso s J 

where we have used (gauge invariant) Schwinger regularisation. The renormalization proce- 
dure we use corresponds to the zero momentum subtraction scheme. Taking the parameter 
So — > generates the ultraviolet divergence. The constant 60 is nothing but the coefficient 
of the beta function to lowest order. Its value arises as a sum of the divergent parts of the 
three terms in Eq.(|2^) which give for Uf" , U^" and Ul;" contributions 3/16, 1/48 and 1/4, 
respectively, so that 

5o = 11/24 

correctly arises. That we get this is another check on the consistency of our formalism. In 
particular, by renormalising in this way we have combined the 0{B^) quantum corrections 
with the classical term, so that the next corrections begin at 0{B'^). 

We thus obtain our final result for the effective potential at high temperature: 

= St) ^ ^WW) ^ °<«"/*^^(^'' + °<»«<^«> ^ P^' 

Since, as T — > 00, M{T) 00 we have gii{M) <ti 1, and our calculation is reliable in this 
regime. So the effective potential acquires a minimum at zero value of the external field. 
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The background field switches off at high temperature, and we can characterise the high 
temperature phase as exhibiting deconfinement. 

4 Discussion 

The central results of this work are expressed in Eqs.( |l3| , ^ ). From these we understand 
that if confinement is due to fields in the QCD vacuum which are long-range (homogeneous 
in our case) and satisfy self-duality or anti-self-duality then this is neither in conflict with 
the Wilson criterion for static quark systems nor with the natural expectation that with 
increasing temperature there is a transition from confinement to deconfinement. It will be 
immediately noticed that we have not obtained a linear heavy quark potential as has been 
observed in lattice simulations. The reason for this discrepancy is straightforward: lattice 
calculations normally implement periodic boundary conditions from the very outset. As 
shall be repeatedly seen in the following, the existence of fields non-vanishing at infinity 
entails significant problems for incorporation on the lattice due to the quasiperiodic bound- 
ary conditions. It seems that there is strong evidence that lattice calculations of the heavy 
quark potential have quite correctly not seen a quadratic potential because the effects of the 
vacuum field we consider have not been built in. How to build these effects in is a problem 
we discuss below. 

On the other hand, the self-dual field, at least at the level of the lowest order propagators 
in this background, does not immediately account for all aspects of confinement: diagonal 
gluons in SU(2) have poles in the propagator and this is a consequence of the fact that 
they do not couple directly to the diagonal background configuration. As mentioned in 
the introduction, the simple self-dual homogeneous configuration is not the entire story, 
and there is room for local effects which can complete the picture of confinement. The 
vacuum field breaks spontaneously CP, colour and 0(3) symmetries. There is a continuum of 
degenerate vacua corresponding to different directions of the vacuum field. This implies ||2^ 
the existence of soliton-like field configurations under the homogeneous background field, 
which could play the role of topologically nontrivial local defects in the QCD vacuum such 
as domain walls. In the absence of explicit solutions we can only speculate on the robustness 
of our results against inclusion of such effects. But insofar as the confining properties of the 
self-dual homogeneous field depend only on the strength of the field and not the direction 
(in real and internal space), it seems plausible that domains distinguished only by changes 
in direction will not disrupt the confinement we observe. 
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However, all of this rests on the assumption that at zero and low temperatures the 
effective potential for this background has a minimum at non-zero field value and there is 
no substitute for a genuine nonperturbative calculation. The only realistic choice for this 
is the formalism of lattice QCD. We thus discuss now in some detail the problems to be 
confronted with setting up the calculation on the lattice, and some insights our preliminary 
investigation into this offers. 

The essential question we need to answer is what the contribution of the homogeneous 
field configuration Eq.(|l]) is to the partition function of lattice SU{2) gauge theory 

Z = J DUeicp{-S[U]} . (35) 
u 

Here, S is now the standard Wilson action and U is shorthand for 

C^n,/x = exp{za^^(an)} e S'C/(2), Vn,/i, (36) 

the link variable, and DU is a functional Haar measure. The lattice spacing is a. Link 
variables are functions of n and are subject to some boundary condition. They thus belong 
to some functional space U. Usually, with representing the size of the lattice in a given 
direction, periodic boundary conditions 

Un+N,fi = Un,fi 

are imposed in order to implement the translation invariance of the theory in the thermo- 
dynamic limit. However, the field B^{an) = aB^yUj^ is evidently not translation-invariant. 
In principle, there are two ways to proceed, both of which have been used in application 
to the Savvidy chromomagnetic background Q. The first choice is to force the long range 
modes to be simply periodic on the lattice. This can be done by 'quantisation' of the field 
strength |, 0, | 

a'B,, = 2vr^, (37) 

where the matrix elements b^y are integers. This certainly provides for periodicity of the 
corresponding link variable, but rewriting Eq.(^7|) as 

B^, = 27r^, L = aN, 
aL 

and going to the thermodynamic (L co) continuum (a 0) limit one obtains 

B^^ = 2-Kb^,y/C, 0<C7<oo, 
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so that the field strength is discretised into multiples of 2tt/C even in the continuum ther- 
modynamic limit. Moreover, the constant C itself depends on our choice of limiting pre- 
scription. These outcomes render this approach rather unappealing. A second approach is 
to change boundary conditions. Free boundary conditions have been advocated in the ap- 



proach of the authors of |1C] who apply it to the lattice calculation of the effective potential 
for the chromomagnetic field in three-dimensional SU{2) theory. 

In our case, Eq.(^) suggests the following generalisation for link variables. We decom- 
pose the general field Afj_ in Eq. (^) into a long range part and the fluctuation Q^: 
A^{an) = Q^{an) + = 1. Quasi-periodic boundary conditions for the 



fields Q can be generalised from Eq. (|18D to all directions now, 

Q^(a(n + N)) = e^"'(")Q^(n)e-™("), w{n) = N^B^pnpT^ ^ /2. (38) 

Thus the following transformation of link variables is generated 

Un+N,t. = e™(") C/n,Me-™^"+'^) • (39) 

This has the structure of a gauge transformation. Thus gauge-invariant quantities such as 
the action are invariant. An integral in Eq. ( |35D includes an integration over all possible 
values of the strength tensor B^y hence all values of w{n), thus the boundary conditions 
Eq.(^) and ( |39[ ) are actually free, consistent with | p!o| . 

The most direct step next is to formulate the effective potential as the lattice functional 
integral 

/ DUe^p{-S[U -V]} 
Juq 

with V denoting a link variable generated by the background field, 



K,^ = exp [ia^B^,nyf/2] 



and where Uq is now the space of quasiperiodic functions. Here B^jj and are external and 
particular directions in the color and euclidean space can be fixed. The actual problem is 
to find an appropriate representation of the measure of the integral such that the exclusion 
of the given background field is manifest. 

The consequence of the gauge function w{n) in Eqs. ( ^J39| ) being nonzero over the 
whole lattice is that all degrees of freedom are affected by the gauge transformation. Thus 
inclusion of covariantly constant field configurations in the space of integration lA in Eq. ( |35|) 
means actually that the space of allowed gauge functions cannot be restricted to the class 



18 



of functions with local support. A significant consequence of this is that Elitzur's theorem 
which forbids spontaneous breakdown of local gauge symmetry [30| does not apply to this 
situation. This theorem concerns the integral 

lim lim {F{U))n,j= lim lim f DU f DgF{U3) exp {-S[U] + JU^} , (40) 

u g 

where F{U) is gauge noninvariant, and J is an external source which breaks gauge invari- 
ance. The order of limits is important. The theorem states that if gauge transformations 
Q are local - namely that they act on a finite (independent of A^) number of degrees of 
freedom - then for sufficiently small sources || J|| < e the following inequality holds 

|exp{JC/5} _ i| < ry(e) , (41) 

with r/(e) being independent of N and vanishing as e goes to zero. Periodicity of the 
functions in li is implicit. 

As has been argued above, both conditions exclude covariantly constant field configura- 
tions which are long range modes that can produce symmetry breaking. Periodic boundary 
conditions and locality of gauge functions are in conflict with a self-consistent incorporation 
of these modes in the lattice functional integral. The choice of free boundary conditions for 
lA and, in particular, the presence in Q of gauge transformations which can act on all degrees 
of freedom results in nonuniformity in the function r/ in lattice size A^, so that JU becomes 
an extensive quantity. In view of this, the drastic difference in the results of (some 
evidence of nontrivial minimum with free boundary conditions and continuous field) and 

1^ (minimum at zero field strength with periodic boundary conditions and 'quantised' 
field strength) seem unsurprising. 

It would be instructive to give an example illustrating that inclusion of the homogeneous 
fields into the lattice integral allows the existence of an order parameter that is not gauge 
invariant. Let us consider the integral over li and Q which includes now homogeneous fields 
and gauge functions of the form Eq. (^8|) . If we put 

F{Ua) = Imi7„3,^„ 

where Un,^u is a plaquette variable, and choose the source term in the form 

^ J^j^yTi Imr^t/^^^^, J^y = const, 
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then the inequahty is not uniform in N for all field and gauge functions: if Un, 



contains the long range fields a B^yUyf t /2 and gauge transformation correspond to 
uj{n) = a^NaB'^pHpf'^T^ /2, then one gets for the source term 

^ J^^Tr inie-*^(")r3e*^(") eY.^{-ia^ B f r^} = 
-2 ^ J^, sin(a2i?^,) [f - 2(/'W' - f ) sm\a^No.B'^f,np/2) 

n,fiu 

Let for simphcity B'^^ = = B'23 = B'^^ = 0, = -^^4 = B' , and Ni = N2 = N'i = 
= N. Using summation formulae 

N 

sin^(nx) = N/2 — cos{N + l)xsin A^x/2sinx, 

n=l 

N 

cos^(nx) = N/2 + cos(A^ + l)x sin Nx/2 sin x, 

n=l 

A , , iV + 1 TVx X 

> sm(nx) = sm x sm cosec— , 

2 2 2' 

n=l 

one gets in the limit — > 00 

AT 

sin^ [a^B'iV(ni - n2 + na - 714) /2] 

ni, 712, ".3,714=0 

iV ^ 

= y J2 {sm'^[a^B'N{n3-n2-ni)/2]+cos^[a^B'N{n3-n2-ni)/2]} + 0{N^) 

712,713,n4=0 

= iV^/2 + 0(iV3). 
Thus for the source term we arrive at the result 

_Ar4 (4/3 _ 2f'^ff'') Y V Ma'B^,) + 0{N'), 

which shows that the gauge dependent part of the source term JU is an extensive quantity, 
and the order of limits J ^ and N ^ 00 cannot be interchanged. 

It should be stressed that this example in no way violates Elitzur's theorem, but just 
underlines that its conditions are too restrictive for a self-consistent incorporation of ho- 
mogeneous field configurations into the lattice functional integral (as mentioned also in the 
last reference of |l30||). 

We repeat that the picture of confinement with a self-dual homogeneous field can become 
reliable only with the inclusion of domain structures in the vacuum such that the symmetries 
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broken by this field are restored at the hadronic level. The boundaries of the domains should 
be describable by some solitonic classical configurations. As far as we are aware, appropriate 
solutions are unknown. It thus remains a problem to verify our considerations of the Wilson 
criterion and Elitzur's theorem in the presence of domains, though we have given plausibility 
arguments why our results might be unaffected. The conclusion is thus that there are two 
interesting unsolved problems to be confronted which may be significant for understanding 
QCD vacuum structure: a calculation of the effective potential for the field Eq.(||) in the 
strong coupling limit, and a search for topologically non-trivial classical configurations in 
the background of a homogeneous self-dual field. 
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Appendix A: Heavy mass limit. 



Below we derive Eq. ( pUj ) starting with Eq. (|6|). In the integrand of Eq. (|6|) let us integrate 
over ^4 

J De4exp|- JdT ^ + ie4(r)(ie(r)) I =e-i^exp|- J dT^iEiir)) 

where T = — 1/4 and we have changed variables 

a = T/ms, T = f3/ms, ^ — ^ 



/ms 



Then the propagator takes the form 

00 

S{x,y\B) oc j dse-'^('+^) j Z)r/exp <j - j 



T 



dp 



2 2ms 2[rnsY 
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where we have omitted the phase factor and a constant in front of the integral. One can 
see that for T ^ oo or more precisely for 

\x — y\ 

and \B\ ~ \E\ <^ m? the integral over s can be evaluated using saddle-point approximation. 
The saddle-point is s = 1, hence we arrive at 



5(f , y, T\B) oc e"*^^ j Dffexp \- J dP 

{ 



^ff-'-B[ffxv\ + ^{Eff? 



(A.l) 



Inserting this result into Eq. ( [A.l| ), we arrive at the representation Eq.(|lO|). 
At zero field E = B = we arrive at the correct nonrelativistic limit 

5(x,y,r|S) oc e-'^^-f ^ =exp|- + r| , (A.2) 

V = « 1, (A.3) 

where the energy of the particle is 

E = m + 



2 

Appendix B: Static Temporal Gauge and Self-Dual Fields 

Consider a fixed direction in colour space such that the phase of the gauge invariant Polyakov 
loop TrP exp(ig /^f dx^Alr" /2) is in the direction. The background self-dual field is also 
chosen to point in the same colour direction. Gauge transformations now used to fix the 
gauge further may not change colour axes. 

An arbitrary gauge transformation on the gauge field = j4^t"/2 takes the form 

4 ^ UA^U^ + -Ud^U^ = A^ . (B.l) 

Under the decomposition of into background B^ and fluctuating parts we choose 
the separate pieces to transform under 'quantum' gauge transformations: 

B^ = UB^U^ 

= UQf^U^ + -Ud^U^ . (B.2) 
9 
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Since we have specified the direction of the background, = 

We use the type Eq.( |B.2D to achieve a fixing of the gauge on Q^. These fields satisfy 



quasiperiodic boundary conditions, Eq.(15), 

Df = 5''^di-e^''^Bi (B.3) 

which must not be changed under gauge-fixing. Thus, U must be quasiperiodic. The 
temporal (axial) gauge (^4 = cannot be achieved with such a group element. The static 
temporal gauge Eq.(|T4|) followed by diagonalisation of the zero mode is however possible. 
Explicitly the U bringing an arbitrary into this gauge is 

[7[Q] = e-*3^4g4r /2pg^p(.^ / dtQi{t,x)) (B.4) 

Jq 

with (?4 = ^ /o^ dtQ/^. Since is quasiperiodic so too is U . Thus are also quasiperiodic: 
e^^*Q^{xi,x) = Q^{x4,x). 

There are still two classes of gauge symmetry remaining: 1. Transformations V{x) = 
exp[z5rt(;^(3;)r^/2] with ^^(x) strictly periodic in x^ can be fixed by introducing an extra 
Lorentz/Coulomb gauge condition on the zero modes of the remaining neutral fields. As 
it is Abelian, this gauge fixing does not introduce Faddeev-Popov ghosts. The one-loop 
effective potential considered in the main body receives no contributions from these neutral 
fields. 2. Transformations M^(x) = exp cause a shift of 2mr/gP in the zero 

mode of the fluctuating field 54. These are relevant for what follows. 

We now implement these considerations in the quantum theory, using the Faddeev- 
Popov trick in the functional integral. The Faddeev-Popov determinant is defined by 

A^'iQ] = J Vg6[Fm (B.5) 

with all configurations related by gauge transformations Eq.(B.2) to a representative 



configuration Q which satisfies F[Q] = 0. The functional F that selects this gauge is 
independent of the background field (unlike in the background field gauge). Inserting 
unity into the generating functional we obtain 

=AfJ VQ^VgS[F[QS]]AF[Q] exp(-5[S + Q]) . (B.6) 



Now we perform a gauge transformation of type Eq.(B^) to bring Q. The measure 



and determinant are invariant, as stated. To recover the same action, a corresponding 
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rotation of the background field must take place, as in Eq.( |B.2D . Because Z is ultimately 
a functional only of the gauge invariant combination S^, we recover again the same Z. We 
may now absorb the integration / Dg into the normalisation in the usual way and obtain 

zyB"] =AfJ VQ^6[F[Q]]Af[Q] eM-S[B + Q]) . (B.7) 

The form of the determinant for the static temporal gauget is well known. Using a lattice 
regularisation for space x, it can be written as 

AF[Q]=llsm\gM^)/2). (B.8) 

X 

The Jacobian is independent of the background component B4. The zeroes of the Jacobian 
indicate the appropriate range of integration for 54, which in turn is seen in the symmetry 
under transformations W at the classical level. The appropriate functional integral after 
implementing the delta functional is then 

VQ,{x) P(?4(x)sin2(5/3g4(f)/2)exp(-5[B + Q]p[Q]=o) . (B.9) 

This is still symmetric under the V transformation. Performing the Faddeev-Popov trick 
again with the Lorentz/Coulomb gauge condition on the neutral zero mode fields enables 
factoring out of this redundant gauge volume. With the normalisation J\f being done at 
-6 = 0, the neutral field contributions to this functional integral will anyway cancel out. 
The W symmetry is however fixed by restriction of the range of integration of ^4. 
We now show how the integral over (74 can be performed. We consider 

j PQiPg4sin2(c//3g4/2)exp|-5 + j (fxJQ^ . (B.IO) 

We integrate over 54 in a diagrammatic expansion in order to derive an effective theory for 
Qi. Thus fields Qi and appear only in external lines of the diagrams. This is a strong 
restriction on the allowed diagrams. The zero mode §4 couples only to charged gluons via the 
three- and four-point vertices, and never to itself. The three-gluon vertex leads to a 54 ^ 
two-charged gluon vertex, while the four-point vertex gives (94)^ — > charged-anticharged 
spatial gluons. This means that the perturbation series in this functional integral stops at 
one loop. Only three topologically distinct classes of diagrams are present and of these only 
one specific diagram gives a non-vanishing contribution as the lattice spacing is taken zero 
after subtraction out of a pure infinite constant. This leaves a mass term in the off-diagonal 
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fields. Its form is determined by the propagator for two fluctuating fields 94 (x), namely 

{0\T{qi{x)q4{y))\0) = . 2. ^ f^.,^. -Sa\B+0] 

J VQi j Vqi{z) sm.'' [gf3qi{z) /2)e ■sol^+^l 

with representing the action with the couplings between 54 and the remaining fields 
dropped. The unregularised form for this is 

So[B + Q] = d^x{Bi{x) + qi{x))V\Bi{x) + qi{x)) 

= d^xqi{x)V^qi{x) 

= So[Q] (B.12) 

because V^B^ = didiB^jXj = diB^i = and q4{x) — > as |x| — > 00. We see that the 
arguments of [|l6| go through unchanged: we shift by half the fundamental domain (74(2?) — > 
qi{x) — tt/qP = q'/^{x) so that the Jacobian becomes a cosine-squared and the boundary 
conditions in X4 of the charged gluons acquire an extra term: they go from being quasi- 
periodic to quasi- aniiperiodic. Next we discretise x = In, with directional unit vectors e, 
and dimensionless field ipfi = g(3q'^{x). We obtain for the action 

Sq\Q] = ^ V>n[^n+2e " '^^n+e + ^n] ■ (B.13) 

Thus the weight factor appearing in the functional integral is 

e-''o'M = f:Cri4^Y . (B.14) 
The functional integrals in Eq.( B.li] ) can be done explicitly 

So, also in the presence of the self-dual background field, the correlator of the fluctuating 
part of Polyakov loops is ultralocal, being proportional to 6^^\x — y) in the continuum limit. 
The result Eq.( p.l^ ) guarantees that the mass term for the charged gluons is as derived by 



(0|r(g4(x)g4(y))|0) = 7-2752 - 2)<5„^^^. + Oi—^) . (B.15) 



ig, namely Eq.(|T5D. 

We are thus led to an effective action after integration out of 54 which contains charged 
1 2 

gluon fields Q^' with a mass diverging with increasing temperature but, in the presence of 
the self-dual field, quasi-antiperiodic boundary conditions. Moreover, the background field 
component i?4(x) is still present in the action in the usual terms where the original A4 was 
located, but the field (74 has been successfully integrated out. 
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Appendix C: Orthogonality and Completeness Relations. 

In this section we derive the orthogonahty and completeness relations for the eigenfunctions 
of the Laplace operator in the presence of the self-dual homogeneous field and the formulae 
for the propagator subject to quasi-antiperiodic boundary conditions. Let us recall first the 
solution to the eigenvalue problem at zero temperature ^ 

-V'^{x)7p{x) = X^{x), 

V\x) = (d^ - iB^{x))\ = ^B^^x, 

in the space of functions vanishing at infinity. The operator — can be represented in the 
form: 



-V2 = 2B (alQ-,a, + 1 



1 /I „ „ \ + 1/1 



a,, = I -Bxu + dn I , a 



Q% = (V ± ib^,) /2, Q±Q± = Q±, Q^Q^ = 0, 6^, = B^,/B. 

The matrix {ib^y) can be diagonalised by means of an appropriate unitary transformation 
U: 

f/^ift;/ = diag(l, -1, 1, -1), C/^a = a, [a^,al] = 5^u, 

ai = (oi + ia2)/\/2, 02 = (ai - ia2)/\/2, 03 = (03 + ia4)/\/2, = {a^ - ia^) / \/2. 



The eigenvalue problem then takes the form: 

2B 1^0202 + 0404 + 1^ il){x) = A^(x), 



with the solution 

\k^k2 = 2B{ki + k2 + l), 

tpkikikski = ^^^-^ (Q2 ) (Q4 ) i^OOks ki{x), 

^ook,k, = ^-j^{a\)'^{al)'^exp [-\bx') ■ (CI) 
The orthogonality and completeness relations have the following form: 

/oo 
d^xi'k,k2k3kS^)^k[k'^k'^k'Sx) = 5k^k[h2k'Jk3k!jkik'^, 
-00 

53 'PLk2ksk4i^)i^kik2k3kAy) = ^ix-y)- (C.2) 

kik2k3k4 
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One sees from Eqs. ( p.lD that the spectrum of is infinitely degenerate which is a conse- 
quence of the homogeneity of the background field. To proceed further, it is advantageous 
to introduce the eigenfunctions (/)fc^fc2(x, y): 



(B /2)'^^~^^* f 1 \ 

(I}kik2{x,y) = ^ E y — — (2/1 - iy2f''{y:i - iyif" exp y--By'^J ipkik2k:ikAx) 

Now the degeneracy is parametrised by the continuous variable y. The function (pQQ{x,y) 
can be seen as a matrix element of the projector onto the subspace spanned by the lowest 
mode {ki = /c2 = 0). 



Making use of Eqs. ( C.2 ), we arrive at the following equations for the case of infinite (3 
(zero temperature): 

-V^(2;)A(x,y) = 5{x - y) 
-V^(a;)(/)fc^fc2(x,y) = Xkik24'kik2{x,y), 
<Pkik2ix,y) = 4>kik2ix - y)exp{-x^5^,,y,,}, 

/oo 
d'^y4>k,k2(^'y)^kik2{z,y) = s{z - x), 

00 

d'^x(t'kik2(^^y)'^k[k'2(^,z) = ^kik[h2k'^4>m{y,z). 

-00 



Together with Eq. (C.3), these define respectively the Green's function and eigenfunctions 
for the operator as well as completeness and orthogonality relations for the eigenfunc- 
tions. The propagator can be decomposed into a sum over projectors onto the subspaces 
corresponding to the different eigen-numbers 

Vk^k2{z,x) 



A(z,x) = 



\2 

kik2 fcife 



Vk,k2 {z, x) = / (fy(l)k,k2 {z, y)4>l,k2 y) = ^fcifc2 {z - x)e5^"^''-"'^ . (C.3) 
J —00 

The completeness, for instance, is derived in the following way: 

/oo 
c^^y E 4>k^k2^x^y)^k^k2{z,y) = 



kik2 

E ^k^k2hikS^)^k^k2k-,kAz) 

kik2kzk4. 
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= Yl '^l^k2kzkS^)'^k^k2kzkAz) = - x). (C.4) 
kik2kzki 

Then at finite /3 tlie function 

CO , • 

D^{x,y)= (-irA(x4-y4 + n/3;x-y)cxpUx^5^^y^ + m/?S4(f + (C.5) 

n=— oo 

is a solution to the equation 

-V\x)D^{x,y) = 5p{x,y), 

Sp{x,y) = Y i-'^TS{x4-yA + nP)S{x-y)expi.-Xij,Bi^,,y^ + inPB4{x + y)\, 

n=— oo 

satisfying the boundary conditions 

-D^(x4 + f3, x; y) = -D^{x4,, x; y) exp {-i/354(f )} , 
D^{x- y4 + P,y) = -Df^{x : y^, y) exp {il3Bi{y)} . 

Here 5(}{x^ y) is the (5-function on the linear space of functions f(}{x) obeying the bound- 
ary condition 

fl3{xA + I3,x) = -//3(x4, x) exp {- 1/3^4 (f)} . 



/oo rfS 
d^y / dyi5p{x,y)fp{y) = fp{x). 
-oo Jo 



One can check that completeness is satisfied: 

-oo ^0 

Moreover, the functions 

^kik2(^^y)= J2 i-^T(pkik2{x4-y4 + nl3,x-y)expl-Xi^Bi^^y^ + inPB4^{x + y)> 

n=— oo ^ 

being the eigenfunctions of (x) , 

and satisfying the boundary condition 

<?^fifc2(^4 + /3,^;y) = -(l)'l^i,^{x4,x)exp{-ipB4{x)} , 



28 



then give 



d^x j dx4(l)^^l^^ix, y)(l)f^i^ix, z) = Sk^i^5k2i^(l)^Q{y, z) 

-oo JO 

SO that they define an orthogonal complete basis for the space $^3. The propagator can 
be decomposed over projectors 



d^z,x) = y: 



\2 ' 

kik2 kik2 

00 13 

^fcifca^^'^) = / J dyA(l)k^k2{z,y)4>i^k2^x,y) 
-00 

= ^ {-l)'^Vkik2{z4- X4 + nP,z- x)ex.p I -Zfj_Bfj,i,Xi, + imj3B4{z + x)> 

n=— 00 ^ ■' 

Taking into account the representation of the zero temperature propagator in terms of the 
projector operators we get Eq. ( p.5| ). 

Let us show how the completeness of the set can be derived. We have to evaluate the 
integral 

-00 r0 ai. o °° r+°° rP °° 



E / dh dy,4l^{x,y)4>l,^{z,y)= ^ / d^y dy, ^ (-1) 
^^k^k2^^'^ - y4 + nP,x - y)4>kik2{z4 -y4 + mP,z- y) 



X exp <y--Xf,Bf,yyy + ^z^B^^yv - infiB^ix + y) + imfiB4{z + y)| . 

After the change of integration variable 1/4 = y4 — nP we get for the right hand side of this 
equation: 

°o r+00 °o /-(l-n)/? 

E / 'i'y i: (-1)"+'"/ dy', 

feifc2=0 °° n,m=-oo -^-^P 

><4ifc2(^4 - y4,x-y)(t)kik2{z4. -y4 + {m- n)P,z- y) 
X exp + ^z^B^^yl + i{m - n)l3B4{z + y)| . 

Finally, shifting the variable m' = m — n in the sum and denoting z'^ = Z4 + m! (3 we arrive 
at 

d'y / 

-00 JO 



E d^y dy4(p1lk^{x,y)(t)^^^^Jz,y) 



felfe2=0 
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°^ ; °° f + OO POO 

= J2 (-l)"" H / '^^y / dy4(l)l^^^{x,y)(l)kik2{z':y)ew{im'PB4{z)} 

m'=-oo feife2=0"^-°° 

= ^ (-l)"'(5(z4 -X4 + m^)(5(z-f)exp<|-2:^S^,,x,, + ?m^B4(z + f)|- = 5^(z,a;). 

m=— oo 

The orthogonality and decomposition of the propagator over the projectors can be obtained 
in a similar way. 

Appendix D: Gluon Propagator in Self-Dual Background Field 
in the Temporal Gauge 

We start with the function 

Dij{x,y\B,m) = Aij{x4 - y4;x - y\B,m)exp 
which is a solution of the equation 



Djk{x,y\B,m) = 6ikS{x - y) 



(- + 5ij + DiDj - 2iBij 

in the limit of infinite /3. The function A is then a solution to the equation 

(-V=^ + dij + DiDj - 2iBij\ l^jk{x\B, m) = 5ik8{x). 

The matrix iBij can be diagonalised by an appropriate unitary transformation U, so that 
we arrive at 



(- + m^) 6rs + DrlD', - 2B5rsis\ ^st{x\B, m) = 8rt6{x) 



(D.l) 



with r, s,t G {0,1,-1} and = s the gluon spin projections onto the third spatial axis. 
Moreover, 

V^ = Dl + D',Dr, D, = uljDj, b', = DjUjs. 



We next decompose the propagator as 



A 



rs — ^rs 



(D.2) 



Using this and the relations 

[Dl Ds] = 2iB6soD4, [Dr, D'.Ds] = 2B5rtitDu 
[D4, Ds] = [D4, D'^] = iB5so, [Do, Ds] = 0, 



^ [-DX, + DrDi - 2B5rs^s 



Ds = 0, 
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we can rewrite Eq. (D.l) as a system of differential equations for the functions H, L, N, 
P: 

(_V2 - 2BCs + m^) F,{x) = 5{x), 

WoLsix) + + m^) Hs{x) + = 0, 

-2iBDiHs{x) + (-V^ - 2Bi, + m^) L,{x) = 0, 

-iD'oPsix) - 2iBD4Hs{x) - + m^) iV,(x) = 0, 

2iBDi {Ls{x) - Ns{x)) + 2B'^Hs{x) + (-V^ + m^) = 0. (D.3) 

We now show that, for the calculation of the space-time trace of the propagator, we only 
need to know the functions H, L, N and P in the neighborhood of x = 0. Consider the 
integral 

r rP 

Tjk= / d^x / dx4D'^u{x,x\B,m), 
Jv JO 

which is contained in Eq. (p^. With L3 the length of the third space direction and using 



Eqs. (I^i^) we get 

Tjk= f d^x 1^ dx4 £ Ajk{n(3■,0\B,m)exp(T^nPBx3^ 



V(3 



Ajk{0;0\B,m) +4 lim V (-l)''Ajfc(n/3; 0|S, m)- 



sin ( ^nBL^P 



n=l 



nBLsP 



VP [Ajfc(0; 0\B, m) + (ij 



(D.4) 



Equation (D^) leads to the result that the terms with n 7^ do not contribute to the 
effective potential. Further calculations can be simplified due to this property of the external 
field Eq.(|D. The solution of Eqs. (^) depends on X , and the first order derivatives are 
proportional to xq = UqjXj or X4. If we need the propagator only for x ^ 0, we can omit 
all terms which contain the first order derivatives. Thus, in the limit x — > 0, we have to 
solve the equations 

(-V^ - 256 + m^) Fsix) = 5{x), 
(^-Dj + m^) Hs{x) + Fs{x) = 0, 
(-V^ - 2B^, + m^) Lsix) = 0, 
(^-Dj + m^) Ns{x) = 0, 
(-V^ + m^) P,(x) + 2B^Hs{x) = 0. 
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One can check that the positive-definiteness of the spectrum of the operators {—D\ + m?) 
and (-V^ - 2Bis + m^) in the space of functions vanishing at infinity. This means that 
Ls{x) — > and Ns{x) for x ^ 0. Finahy one gets for small 

Fs{x) = (^-V^-2B^s + m^y^5{x), 
Hs{x) = -(-Dl + n?y^ Fs{x), 
Ps (x) = -2^2 (- + m^) Hs (x) , 

Ls{x) = 0, Ns{x) = 0. (D.5) 

Using Eqs. ( p.5| ) and ( p.2| ), we arrive at the relation Eq.(p9[) given in the main body of 
the paper. We have only to insert the delta-function into the equations for Hs and Pg and 
to represent these functions as convolutions of the propagators F^, Fq and A4, where the 
latter is the Green's function corresponding to (— -D4 + m^). These lead to the expressions 
in Eq.d^). 
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